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We  analyze  the  minimum  achievable  mean-square  error  in  frequency-modulated  continuous-wave  remge 
estimation  of  a  single  stationary  target  when  photon-counting  detectors  are  employed.  Starting  from  the 
probability  density  function  for  the  photon-arrival  times  in  photodetectors  with  subunity  quantum  effi¬ 
ciency,  dark  counts,  and  dead  time,  we  derive  the  Cramer-Rao  bound  and  highlight  three  important 
asymptotic  regimes.  We  then  derive  the  maximum-likelihood  (ML)  estimator  for  arbitrary  frequency 
modulation.  Simulation  of  the  ML  estimator  shows  that  its  performance  approaches  the  standard  quan¬ 
tum  limit  only  when  the  mean  received  photons  are  between  two  thresholds.  We  provide  anal3dic  approx¬ 
imations  to  these  thresholds  for  linear  frequency  modulation.  We  also  compare  the  ML  estimator’s 
performance  to  conventional  Fourier  transform  (FT)  frequency  estimation,  showing  that  they  are  equiv¬ 
alent  if  the  reference  arm  is  much  stronger  than  the  target  return,  but  that  when  the  reference  field  is 
weeik  the  FT  estimator  is  suboptimal  by  approximately  a  factor  of  \/2  in  root-mean-square  error.  Finally, 
we  report  on  a  proof-of-concept  experiment  in  which  the  ML  estimator  achieves  this  theoretically 
predicted  improvement  over  the  FT  estimator.  ©  2013  Optical  Society  of  America 
OCIS  codes:  010.3640,  270.5290,  030.5260,  000.5490,  040.1345. 


1.  Introduction 

In  its  essence,  optical  ranging  is  a  problem  of  esti¬ 
mating  the  round-trip  flight  time  of  a  phase-  or 
amplitude-modulated  optical  beam  that  reflects  off 
of  a  target.  Frequency-modulated  continuous-wave 
(FMCW)  ranging  systems  obtain  this  estimate 
by  performing  an  interferometric  measurement 
between  a  local  frequency-modulated  laser  beam 
and  a  delayed  copy  returning  from  the  target.  The 
range  estimate  is  formed  by  mixing  the  target-return 
field  with  the  local  reference  field  on  a  beam  splitter 
and  detecting  the  resultant  beat  modulation.  In 
conventional  FMCW  ranging,  the  source  modulation 
is  linear  in  instantaneous  frequency,  the  reference- 
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arm  field  has  many  more  photons  than  the  target- 
return  field,  and  the  time-of-fiight  estimate  is  gener¬ 
ated  by  balanced  difference  detection  of  the  beam 
splitter  output,  followed  by  a  frequency-domain  peak 
search.  This  technique  is  known  to  achieve  range  res¬ 
olution  proportional  to  c/2A/‘  in  a  T-second  chirp 
window,  where  c  is  the  speed  of  light,  Af  is  the  chirp 
bandwidth,  and  the  time-bandwidth  product  satis¬ 
fies  TAf  »  1  [1,2]. 

Important  advances  have  been  reported  in  the  re¬ 
cent  literature  for  FMCW  ranging  systems.  To  revisit 
several  examples,  very  broad  chirp  bandwidths  have 
been  achieved  using  active  feedback  stabilization 
[3,4].  Alternative  frequency -modulation  functions, 
such  as  sinusoidal  modulation  [5,6]  and  pseudoran¬ 
dom  modulation  [7] ,  have  been  employed  to  improve 
the  range  resolution.  These  results,  and  much  of 
the  literature  in  general,  have  focused  on  the 
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strong-local-oscillator  regime,  wherein  the  mean 
photon  number  in  the  reference  field  greatly  exceeds 
that  in  the  target-return  field  over  the  integration 
window  of  interest.  A  strong  local  oscillator  provides 
optical  gain  of  a  weak  target  return  prior  to  detec¬ 
tion,  such  that  the  predominant  noise  in  the  photo¬ 
current  is  the  quantum  (vacuum)  noise  of  the 
incident  field  mode  [8,9]. 

With  high  efficiency,  high  bandwidth,  and  low- 
noise  photon-counting  detector  technologies  rapidly 
maturing,  it  is  feasible  to  anticipate  that  FMCW 
ranging  systems  could  utilize  photon-counting  detec¬ 
tors.  These  detectors  offer  very  high  sensitivity,  and 
can  accommodate  the  large  bandwidths  needed  for 
ranging  applications,  but  their  limited  dynamic 
range  (i.e.,  saturation)  and  their  susceptibility  to 
dead  time  (i.e.,  blocking)  typically  precludes  the 
use  of  a  very  strong  local  oscillator.  In  recent  years, 
there  has  been  interest  to  quantify  the  performance 
attainable  with  high-sensitivity  photon-counting 
detectors  paired  with  a  local  oscillator  of  strength 
comparable  to  the  target-return  signal  [W] — that 
is,  a  weak  local  oscillator.  In  this  paper  we  continue 
to  investigate  this  question,  focusing  both  on  deriv¬ 
ing  an  asymptotically  tight  lower  bound  to  the  esti¬ 
mation  accuracy,  and  on  identifying  the  performance 
of  an  estimation  algorithm  that  approaches  this 
bound.  Because  range  is  related  to  time  delay  by  a 
constant  (the  speed  of  light),  we  present  all  of  our 
results  in  terms  of  the  latter  and  omit  the  trivial  scal¬ 
ing  of  these  results  to  express  them  as  a  range.  Our 
treatment  is  founded  on  a  rigorous  statistical  charac¬ 
terization  of  the  (random)  photoelectron  emission 
times  as  a  function  of  the  incident  optical  field, 
including  the  deleterious  effects  caused  by  dark  cur¬ 
rent  and  dead  time.  These  statistics  permit  us  to  de¬ 
rive  the  Cramer-Rao  lower  bound  (CRB)  on  the 
accuracy  of  FMCW  ranging  and  to  derive  the 
maximum-likelihood  (ML)  estimator,  whose  perfor¬ 
mance  approaches  this  bound  under  some  photon- 
flux  conditions.  Here,  we  determine  these  conditions 
and  validate  the  theory  via  simulation  and  experi¬ 
ment.  Because  our  emphasis  is  on  FMCW  receivers 
employing  photon-counting  detectors,  we  limit  the  fo¬ 
cus  of  this  paper  to  analyzing  the  performance  when 
a  single  specularly  reflecting  target  is  present  within 
the  range  of  interest.  Nonetheless,  the  formalism  de¬ 
veloped  herein  could  be  extended  to  include  diffuse 
target  reflections  [11],  multiple  targets  within  the 
range  of  interest  [12],  and  atmosphere-induced 
distortions  [13],  with  the  appropriate  modifications 
to  statistics  governing  the  output  processes  of  the 
photon-counting  detectors. 

Our  paper  is  organized  as  follows.  In  Section  2  we 
define  the  problem  of  interest,  and  then  we  derive  the 
CRB  for  arbitrary  frequency  modulation.  First,  we 
consider  only  the  shot  noise  of  the  incident  (coherent- 
state)  fields  and  determine  the  CRBs  achieved  by  dif¬ 
ferent  modulation  functions.  We  then  extend  our 
formalism  to  include  dark  counts  and  dead  time 
in  the  detector,  and  quantify  the  degradation  in 


performance.  In  Section  3,  we  derive  the  ML  estima¬ 
tor  and  show  via  simulation  that  when  the  total 
mean  received  photon  number  is  between  two 
thresholds,  its  mean-square  error  (MSE)  approaches 
the  standard  quantum  limit.  We  derive  analytic  ap¬ 
proximations  to  these  threshold  values  for  the  case  of 
linear  frequency  modulation.  In  addition,  we  com¬ 
pare  the  ML  estimator’s  performance  to  conventional 
frequency-domain  estimation,  highlighting  regimes 
in  which  they  perform  identically  and  regimes  in 
which  they  do  not.  Finally,  in  Section  4,  we  report 
on  an  experiment  with  linear  frequency  modulation 
and  compare  the  results  to  our  analysis. 


2.  Problem  Formulation  and  the  Cramer-Rao  Bound 

Figure  1  illustrates  the  FMCW  ranging  system  in 
consideration.  We  assume  that  an  ideal  laser  field 
[14]  with  center  frequency  W(.  is  split  via  a  beam  split¬ 
ter  into  a  signal  field  and  a  reference  field 

where  =  ^'N  for 

-T/2  <t  <  TI2,  and  periodically  repeating  other¬ 
wise.  The  subscript  k  =  S,R  denotes  their  complex 
baseband  temporal  modulations,  in  which  the  in¬ 
stantaneous  frequency  modulation  satisfies 
=  0,  where  (pit)  =  d(p{t)/dt.  In  addition, 
Ef^iit)  are  normalized  to  have  units  ^photons/s, 
such  that  Nk  \  denotes  the  mean  photon  number  over 
the  T-second  window.  The  signal  field  is  transmitted 
to  interact  with  the  target,  yielding  a  target-return 
field  Es2{t)-  In  this  paper  we  have  assumed  a 
single  specularly  reflecting  target,  so  Es2it)  = 
rgEs  iit  -  to),  where  <  1  is  the  field  reflection  coef¬ 
ficient  and  ^0  s  ^  is  the  relative  time  delay  between 
the  signal-  and  reference-arm  fields  within  the  a  pri¬ 
ori  domain  of  uncertainty  T.  The  reference-arm  field 
remains  local  to  the  ranging  instrument,  and  is 
therefore  assumed  to  go  through  a  simple  attenua¬ 
tion  by  rR<X  yielding  E^,2(0  =  i(^).  The 

receiver  mixes  the  target-return  field  and  the  refer¬ 
ence  field  with  a  50/50  beam  splitter,  resulting  in 


ED,m{t) 


Es.2{t)  +  i-'E)”^^^Eji2it) 
42 
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Fig.  1.  (Color  online)  Block  diagram  of  FMCW  ranging.  Optical 
fields  are  denoted  with  thick  (orange)  lines  and  electrical  signals 
are  denoted  with  thin  (black)  lines. 
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at  its  two  output  ports  m  =  1,2.  are  incident 

on  two  matched  photon  counters  for  t  e  {-T 12,  T /2\ 
giving  rise  to  the  photocurrent  processes  imit).  We 
assume  im{t)  are  normalized  by  the  electron  charge 
to  have  units  photoelectrons/s,  and  are  therefore 
sequences  of  unit-area  impulses  occurring  at  each 
photoelectron  emission  epoch  [15].  With  ideal  laser 
light  incident  on  the  photodetectors,  im{t)  are  inde¬ 
pendent  discrete  stochastic  point  processes  whose 
statistics  are  governed  by  several  uncertainties. 
First,  energy  measurements  on  coherent-state  fields 
are  stochastic  due  to  the  quantum  nature  of  light 
(shot  noise)  [16].  Second,  there  are  stochastic  elec¬ 
tron  emissions  from  the  photosensitive  area  of  the 
detector  when  no  light  is  incident  (dark  counts). 
Third,  after  every  photoelectron  emission  there  is 
a  period  during  which  the  detector  is  unable  to  gen¬ 
erate  any  additional  photoelectrons  regardless  of  the 
incident  field  amplitude  (dead  time).  The  first  of 
these  effects  is  fundamental,  whereas  the  latter 
two  are  practical  limitations  to  state-of-the-art  pho¬ 
ton  counters  that  can  be  improved  upon  in  the  future. 

An  important  component  in  the  Fig.  1  system  is 
the  algorithm  that  generates  an  estimate  tg  of  tho 
time  delay  to,  based  on  the  occurrence  times  of  the 
photoelectrons  from  the  two  detectors.  Different 
estimators  will  achieve  different  MSE  values,  an  ex¬ 
haustive  characterization  of  which  is  not  feasible. 
The  CRB,  on  the  other  hand,  is  a  lower  bound  on 
the  MSE  achievable  by  any  unbiased  estimator  [2]. 
Furthermore,  the  CRB  is  achieved  by  the  ML  estima¬ 
tor  if  an  efficient  estimator  exists,  and  even  if  one 
does  not  exist  the  ML  estimator  often  asymptotically 
approaches  the  CRB  at  high  signal  flux.  Therefore, 
we  first  derive  the  CRB  for  the  FMCW  ranging 
system  described  above.  Then,  we  focus  on  the  ML 
estimator  and  evaluate  its  performance  relative  to 
the  CRB. 

A.  Standard  Quantum  Limit 

Let  us  first  consider  near-ideal  photon  counters  and 
develop  some  insights  into  the  ultimate  performance 
of  FMCW  ranging.  When  and  Ejj2{t)  are 

incident  on  matched  photon-counting  detectors  with 
quantum  efficiency  rj,  zero  dark  current,  and  zero 
dead  time,  im{t)  are  statistically  independent 
inhomogeneous  Poisson  point  processes  with  rate 
functions  [9,17,18] 

cos(D(t;^o)  -I-  lyo)],  (2) 

for  m  =  1, 2.  Here,  D{t\ t^)  =  (p{t)  -  (fft  -  tg)  is  the  beat 
modulation  that  carries  information  on  to,  and  i/zq  = 
Wcto  denotes  the  phase  shift  in  the  carrier.  The 
parameter 

2^Ns,2Nr,2 
'  Ns.2+Nr,2 

is  between  0  and  1,  indicating  the  relative  strength  of 
the  target-return  and  reference  fields  in  terms  of 


their  mean  photon  numbers  Nk.2  =  for 

k  =  S,R,  and  Ni  =  Ns.2+Nr2  is  the  total  mean 
number  of  photons  incident  on  the  two  detectors. 
The  probability  density  function  for  observing  ni 
photodetection  events  in  the  first  detector  at  event 
times  {«!,.. e(-T/2,T'/2]  and  n2  photodetec¬ 
tion  events  in  the  second  detector  at  event  times 
{vi,...,VnJe  (-T/2,r/2],  is  given  by  [W] 


p{ni,Ui,...,Un^,n2,Vi,...,Vn^-,to) 
r  »i  1  r  »2 


~  X^/2  + ^2  (f;fo)] 


(4) 


From  Eq.  (4),  the  CRB  for  to  is  expressed  as  [2] 
CRB  =  1^  iog p{ni,Ui,...,Un^,n2,Vl,...,Vnf,tQ)Y^ 

\  <  I  ’ 

(5) 

where  the  angled  brackets  denote  the  expectation 
over  all  of  the  random  variables. 

Throughout  this  paper  we  assume  that  the  a  priori 
uncertainty  in  the  target  range  is  much  greater  than 
an  optical  wavelength,  such  that  the  carrier  phase 
shift  lyo  provides  negligible  information  in  determin¬ 
ing  the  true  range.  Consequently,  to  obtain  relevant 
lower  bounds  on  attainable  MSE,  we  ignore  the  to 
dependence  of  lyo  by  treating  it  as  a  constant.  In  par¬ 
ticular,  the  theoretical  derivations  below  assume 
lyo  =  0  for  simplicity.  However,  in  applications  that 
have  range  uncertainty  on  the  order  of  a  few  wave¬ 
lengths,  lyo’s  ^0  dependence  should  be  taken  into 
account  when  calculating  Eq.  (5). 

In  Appendix  A  we  show  that  substituting  Eqs.  (2) 
and  (4)  into  Eq.  (5)  leads  to 


CRB: 


1 

rjNi 


If 

tJ- 


■TI2_^^  pHin\D{t-,to))  ^2., 
.T/2  l-p^cos‘^iD(t;to)) 


1-1 


(6) 


Eor  a  fixed  ratio  between  the  target-return  and  the 
reference-arm  mean  photon  numbers,  i.e.,  for  f  con¬ 
stant,  the  CRB  follows  a  l/Nj  scaling,  which  is  re¬ 
ferred  to  as  the  standard  quantum  limit  because  it 
is  the  best  possible  scaling  achievable  with  photon 
counting  of  ideal  laser  fields.  In  addition,  because 
the  first  term  in  the  integrand  is  maximized  at  f  = 
1  for  all  t,  we  find 


CRB  > 


j^n  r 

riNjlTl 


T/2 

T/2 


dt^^it  -  to) 


(7) 


for  all  Nj  >0  and  ^  with  equality  if  and  only  if 
P  =  1.  Several  insights  are  noteworthy  here.  Eirst, 
given  a  total  mean  received  photon  number  Nj,  it 
is  optimal  to  have  P  =  1,  i.e.,  Ns.2  =  Nr2-  Second, 


2010 


APPLIED  OPTICS  /  Vol.  52,  No.  10  /  1  April  2013 


Table  1.  L2  Norms  of  Various  Frequency  Modulations  with  Identical 
Frequency  Detuning  Range 


Type 

ipit),  t  e  (-T/2.  T/2] 

Linear 

Aft/T 

Ap/Vl 

Polynomial 

(Af/2)sga(t)(2\t\/T)P,P>0 

(Aff/(i(2p  +  1)) 

Sinusoid 

(Af/2)  sin(2ffi/r) 

Af^/8 

(Pseudo) 

(Af /2)Yff—K<^i‘bh(t),  Ok  e 

AfVd 

random 

1-1.1} 

“sgn(f)  denotes  the  sign  of  f,  bi^{t)  =  1  for  2|f-  /cA|/A  <  1  and  0 
otherwise,  where  T/A  =  2K  +  1 ,  and  a*  are  independent  euid 
equally  likely  to  take  on  either  value. 


add  dark  counts  and  dead  time  to  quantify  the 
degradation  in  the  CRB  relative  to  the  standard 
quantum  limit.  Suppose  that  the  photon-counting 
detectors  from  Fig.  1  have  equal  and  constant  dark 
count  rates  but  no  dead  time.  Then,  the  dark 
counts  in  each  detector  form  a  homogeneous  Poisson 
process  with  statistics  equivalent  to  having  on  aver¬ 
age  Nd  =  AdT/ri  additional  photons  incident  in  a 
T-second  interval.  Therefore,  the  two  photodetector 
outputs  are  independent  inhomogeneous  Poisson 
processes  as  stated  in  the  previous  subsection,  but 
now  with  the  modified  rate  functions  [21] 


Eq.  (7)  is  independent  of  to  when  the  modulation  is 
periodic  with  T.  Finally,  the  optimal  CRB  depends 
only  on  the  average  energy  (£2  norm)  of  the  baseband 
(zero-average)  frequency-modulation  function. 
Therefore,  frequency  modulations  with  equal  £2 
norm  attain  the  same  optimized  CRB.  Table  1 
summarizes  the  £2  norms  of  several  frequency 
modulations,  quantifying  their  relative  performance 
at  equal  Nj. 

Setting  /i  ^  1  in  Eq.  (6),  we  derive  the  CRB  in  the 
strong-reference-arm  regime  as 


CRB 


1 

WWi 


d^  sin^ {D{t;to))^^ (t  -  to) 


,  (8) 


where  the  term  outside  the  brackets  is  approxi¬ 
mately  equal  to  1/(4;7A^S  2)-  We  can  now  compare 
the  CRBs  for  ft  =  1  and  ft  in  several  relevant 
scenarios.  liNs  2  is  fixed  and  the  reference-arm  pho¬ 
ton  number  is  allowed  to  vary,  the  ft  =  1  and  ft 
cases  yield  approximately  equal  CRBs  given  by 
l/(2riNs.2e,p)  [^1,  where  The 

former  uses  fewer  photons  to  achieve  this  sensitivity 
but  requires  photon-counting  detectors,  whereas 
the  latter  uses  more  photons  but  could  utilize 
linear-mode  detectors.  Another  relevant  scenario  is 
one  in  which  the  mean  number  of  source-generated 
photons  is  fixed,  i.e.,  A^s.i  +  Nr.i  =  W,  and  these  pho¬ 
tons  are  allocated  into  the  signal  and  reference  arms. 
Assuming  the  common  scenario  of  high  loss  along  the 
signal  path  (r|  ^  1)  but  low  loss  along  the  reference 
path  (r|  «  1),  the  CRB  for  the  ft  =  1  case  is  approx¬ 
imately  given  by  l/{2rir\Ne^),  whereas  the  CRB  for 
the  ft  case  is  approximately  y / i2rir'^N e^),  with 

y=l+  Nrj/Ns,i  »  1  -I-  4r|  (9) 

being  strictly  greater  than  1.  So,  allocating  photons 
such  that  ^  ^  1  in  general  loses  relative  to  having 
ft  =  1  for  N  fixed,  but  the  performance  loss  is 
negligible  if  the  mean-photon-number  allocation 
satisfies  1  ^  Nr  i/Ns,i  »  4r|. 

B.  Dark  Noise  and  Dead  Time 

The  results  above  define  the  best  possible  perfor¬ 
mance  attainable  with  ideal  laser  light  and  near¬ 
ideal  photon  counting,  with  subunity  quantum 
efficiency  being  the  only  nonideal  effect.  Now  we 


=  ^[1  +  Cosm-Jo))]  (10) 

for  m  =  1, 2,  where  Njd  =^7-1-  2Nd  and  =  /i/(l  -I- 
2NdlNi)  are  the  two  parameters  that  are  modified 
due  to  dark  counts.  Next,  suppose  that  the  detectors 
also  have  equal  and  deterministic  nonparalyzing 
dead  time  t  [17] .  In  this  case  the  rate  function  itself 
becomes  a  stochastic  process  that  causally  depends 
on  prior  arrival  times.  Consequently,  the  detector 
output  is  modeled  as  a  self-exciting  point  process. 
The  probability  density  function  for  observing  ni 
photodetection  events  in  the  first  detector  at  event 
times  {«!,.. e  (-T'/2,T'/2]  and  n2  photodetec¬ 
tion  events  in  the  second  detector  at  event  times 
{vi,...,VnJ  e  {-Tf2,Tf2]  is  now  given  by  [17] 

pini,Ui,  ...,Un^,n2,Vi,  ...,Vn^-,to) 

=  Pd{ni.Ui,  ...,Un^-,to)Pd{n2,Vn...,Vn^\to),  (11) 

where 

r 

Pdfb'm  1  •  •  •  5  ~  1^1  I  ^o) 

xexpi  -  /  '  dtA^jit-Jo) 

[  J-T/2 

-  /  dtA^4it;  to)wAt  -  w;_i)  \  (12) 

i_2  Jui_i  ) 

for  m  =  1,2.  Here,  we  have  defined  =  T/2  for 
compactness  of  the  expression,  and  wAt)  =  1  if  ^  > 
T  but  is  0  otherwise. 

Deriving  the  CRB  is  straightforward  by  substitut¬ 
ing  Eqs.  (H)  and  (12)  into  Eq.  (5),  so  we  postpone  the 
tedious  details  to  Appendix  A.  Here  we  highlight  sev¬ 
eral  asymptotes  to  the  CRB  expression  when  the  rate 
of  change  of  A^^dilAo)  is  much  slower  than  1/t, 
specifically,  when  the  frequency  modulation  and 
dead  time  satisfy  |^(t)  -  (^(t  -  to)|r  ^  1  for  all  t  El 
{-T/2,  T/2]  and  to  e  T: 


CRB 


Nd 

Afj 

^K'^^ijtAo)  Nj 


«2iVd 

and 

a  1 

»2Nd 

and 

a  1 . 

»2Nd 

and 

a  »  1 

(13) 
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where  a  =  t]Ni  d'^/(2T)  is  the  average  fraction  of  pho¬ 
toelectron  counts  lost  during  the  dead  times,  and  the 
three  terms  independent  of  photon  number  are 


1  pT/2 

Kiito)  =  -p=  (14) 

r  J-T/2 

T  J-T/2  1  -  ^  cos^(D(^;  ^o)) 


(p  (t-to), 
(15) 


and 

Ks{p,tQ)  =  -l  dt 
-t  J-T/2 

sin2(D(^;  ^o))[l  +  cos^{D{t\ ^o))] 

[1- p‘^cos^(D(t\ta))f 

(16) 

The  asymptotic  expressions  in  Eq.  (13)  represent 
three  distinct  trends  with  respect  to  the  mean  inci¬ 
dent  photon  number.  When  the  mean  dark  count  is 
significantly  greater  than  the  mean  signal  photon 
count  {Nfi  »  Ni)  and  the  dead  time  is  short  enough 
to  have  negligible  impact  on  the  arrival  rate  (a  «  1), 
the  performance  is  dark-noise  dominated.  In  this 
regime  the  CRB  is  significantly  worse  than  the  stan¬ 
dard  quantum  limit,  and  it  improves  quadratically  as 
the  mean  incident  photon  number  increases.  At  the 
other  extreme,  when  the  mean  incident  photon  num¬ 
ber  is  so  large  that  the  detector  is  almost  universally 
blocked  by  dead  time  (a »  1),  the  output  of  the 
photon-counting  detectors  saturates  irrespective  of 
the  mean  received  photon  number.  Here  the  perfor¬ 
mance  is  dead-time  dominated  and  the  CRB  reaches 
a  floor.  Often  there  is  an  intermediate  regime  be¬ 
tween  these  two  extremes  wherein  the  signal  counts 
dominate  over  dark  counts  (W/  »  W^),  yet  the  dead 
time  impact  is  negligible  (a  ^  1).  In  this  regime 
the  performance  is  shot-noise  dominated  and  the 
CRB  approaches  the  standard  quantum  limit.  If  dark 
counts  or  dead  time  are  too  high,  the  ranging  perfor¬ 
mance  may  never  become  shot-noise  dominated,  thus 
the  standard  quantum  limit  may  not  be  attainable. 

In  Fig.  2  we  have  plotted  the  CRB,  normalized  by 
to  render  it  dimensionless,  when  j){t)  =  ^ftjT, 
i.e.,  when  linear  frequency  modulation  is  employed. 
The  dashed  curves  correspond  to  the  aforementioned 
asymptotes,  and  the  solid  line  is  the  numerically 
evaluated  CRB  from  Eq.  (A6)  in  Appendix  A.  The 
figure  indicates  that  the  asymptotes  provide  a  good 
approximation  to  the  CRB  in  their  applicable  regime. 
Furthermore,  the  transition  thresholds  are  well- 
approximated  by  the  intersection  points  of  these 
asymptotes. 


3.  Maximum-Likelihood  Estimation  of  Range 

In  this  section  we  derive  the  ML  estimator  and  show 
via  simulation  that  its  performance  approaches  the 


Fig.  2.  (Color  online)  Exact  VCRB,  normalized  by  T,  is  shown  as 
the  solid  (blue)  line.  The  dashed  (red)  line  is  the  dark  noise  asymp¬ 
tote,  the  dash-dotted  (green)  line  indicates  the  shot-noise  asymp¬ 
tote,  and  the  dotted  (cyan)  line  is  the  dead-time  asymptote. 


CRB  when  the  incident  photon  number  satisfies 
<  Nj  <  Nu,  where  Nf  and  are  lower  and 
upper  breakdown  thresholds,  respectively.  The  lower 
threshold  is  due  to  having  a  nonlinear  relation  to 
the  observation  statistics,  whereas  the  upper  thresh¬ 
old  is  due  to  dead  time  dominating  the  observation 
statistics.  We  will  derive  analytic  approximations 
for  Nf  and  in  the  case  of  linear  frequency 
modulation. 

The  ML  estimator  is  defined  as 


to  =  argmax  logp(ni,wi,  ...,Vn^\to), 

(17) 

where  T  is  the  support  of  to-  Using  the  probability 
density  function  in  Eq.  (11)  and  following  straightfor¬ 
ward  simplifications  (see  Appendix  A  for  details 
omitted  here),  we  arrive  at 


/T/2 

dt[ii{t)hx{t-,to)  -\-  i2{t)h2{t-,to)], 

T/2 

(18) 


where  the  impulse  responses  are  given  by 
hmitdo)  =  l0g(l  -h  (-1)'”+Vd  COSiPitdo))) 

+  {-ir+^aPd  cos(D(Lto)).  (19) 

for  m  =  1,  2.  Thus,  a  sufficient  statistic  for  the  ML 
estimator  is  the  sum  of  filtered  photodetector  out¬ 
puts,  where  the  filter  impulse  responses  are  linear 
and  time  varying. 

Let  us  consider  several  limiting  cases  of  this  gen¬ 
eral  estimator.  When  dead  time  is  negligible,  i.e. 
a  1,  only  the  first  term  in  Eq.  (W)  survives,  and 
the  impulse  responses  are  recognizable  as  the  log- 
matched  filters  to  the  photoelectron  generation  rates 
in  each  detector.  If  ^  1  and  a  is  arbitrary,  which 
occurs  with  a  strong  local  field  (/i  ^  1)  or  when  dark 
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noise  dominates  the  signal  energy  (Na  »  Nj),  the  ML 
estimator  can  be  simplified  as 

/TI2 

dt(ii(t)  -  ^2(0)  cos(D(t;  to));  (20) 

T/2 

i.e.,  difference  detection  followed  by  a  matched  filter 
to  the  rate  function  generates  a  sufficient  statistic.  If 
we  further  specialize  our  result  to  ^(t)  =  we 

obtain 

to  «  arg  max 

X  3ft|e''^/''o/<2r)  duihiu)  -  t2(u))e-'^^*»“/^|; 


i.e.,  difference  detection  followed  by  a  Fourier  trans¬ 
form  (FT)  and  peak  detection  (henceforth  referred  to 
as  the  FT  estimator)  is  also  approximately  the  ML 
estimator  in  this  regime.  This  result  is  not  surpris¬ 
ing.  Both  with  a  strong  reference  field  or  with 
dark-count-dominated  noise,  the  photodetector  out¬ 
put  processes  have  a  signal-dependent  mean  but  a 
largely  signal-independent  noise  process.  This,  to¬ 
gether  with  the  (asymptotically  tight)  approximation 
that  the  photodetector  outputs  form  a  Gaussian 
random  process  when  the  photoelectron  counts  are 
large,  yields  an  additive  Gaussian  noise  channel  for 
which  the  FT  estimator  is  exactly  ML  achieving  [2]. 

The  equivalence  of  the  two  estimators,  however, 
does  not  extend  to  the  ^  »  1  regime.  To  show  this, 
we  approximate  the  term  inside  the  brackets  in 
Eq.  (^)  as  S{Q.-  Gp)  -I-  N{Q.),  where  Q  =  Afto/T  is 
the  normalized  variable  over  which  the  maximiza¬ 
tion  is  performed,  Gq  is  the  true  value  of  the 
parameter  with  the  same  normalization,  S(G)  = 
r]Njfi  sin{Q.T/2)/(Q.T)  is  the  mean,  and  N(Q.)  is  a 
zero-mean  complex  Gaussian  random  process  with 
the  nonzero  correlation  function  {N* (^i)N {^2))  = 
{N[n:/T)5{Q.i- 0.2).  Using  the  CRB  for  estimating 
Go  as  the  true  performance  of  the  FT  estimator,  we 
arrive  at  [2] 

((^0  -  «  2CRB  (22) 

for  P  =  X  where  denotes  the  FT  estimator,  and 
CRB  refers  to  Eq.  (6).  Therefore,  as  the  ML  estima¬ 
tor’s  performance  approaches  the  standard  quantum 
limit,  we  expect  that  the  FT  estimator  will  attain  ap¬ 
proximately  a  factor  of  a/2  larger  root-mean-square 
error  (RMSE).  On  the  other  hand,  when  ^  1,  we 
obtain  ((to  -  4^^*)^)  ~  CRB,  where  CRB  is  now  given 
by  Eq.  (8),  in  agreement  with  our  earlier  conclusion 
that  the  FT  and  ML  estimators  are  equivalent  in  this 
regime. 

Figures  3  and  4  plot  the  simulated  MSEs  of  the  ML 
estimator  and  the  FT  estimator  for  linear  frequency 
modulation  as  a  function  of  rjNi.  Figure  3  represents 
the  near-ideal  case  with  =  0  and  r  =  0  and  «  1. 
We  observe  that  the  ML  estimator  tracks  the  CRB 


Fig.  3.  (Color  online)  RMSE  normalized  to  the  width  of  the  range 
uncertainty  window.  The  simulated  ML  (magenta)  and  FT  (cyan) 
estimator  performances,  the  CRB  (dashed  red),  the  anal3dic  RMSE 
model  (dash-dotted  green),  and  the  FT  RMSE  asymptote  (dotted 
hlue)  are  plotted.  Sr  is  the  fractional  position  of  the  object  in  its 
support. 

well  for  large  values  of  Nj.  However,  when  Nj  falls 
below  a  threshold,  which  we  denote  with  N f,  its  per¬ 
formance  rapidly  deteriorates  toward  a  uniformly 
distributed  random  guess  of  to  £  2'.  The  FT  estima¬ 
tor’s  performance  is  similar.  However,  consistent 
with  the  derivation  in  the  previous  paragraph,  its 
RMSE  is  approximately  a  factor  of  \/2  worse  than 


(a) 


riN,  [pe] 

(b) 

Fig.  4.  (Color  online)  RMSE  normalized  to  the  width  of  the  range 
uncertainty  window.  Plot  legend  is  identical  to  that  of  Fig.  3. 
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the  ML  estimator’s  performance  when  Nj  >  Nf. 
Figure  4(a)  includes  dark  counts  and  dead  time,  spe¬ 
cifically  Nil  =  X  T  =  100  ns,  and  ^  «  1.  The  qualita¬ 
tive  behavior  is  the  same  as  before:  for  Nj  <  Nf, 
the  ML  estimator  rapidly  deteriorates  to  random 
uniform  guesses,  and  above  this  threshold  its  MSE 
tracks  the  CRB.  However,  unlike  the  previous  case, 
the  CRB  itself  saturates  when  Nj  exceeds  a  thresh¬ 
old  we  denote  as  Nu,  causing  the  MSE  to  deviate 
from  the  standard-quantum-limited  scaling  of 
l/irjNj).  Between  these  thresholds,  where  the  ML 
estimator  performs  at  the  standard  quantum  limit, 
the  FT  estimator  exhibits  the  V2  gap  in  RMSE,  as 
expected.  Finally,  Fig.  4(b)  shows  a  case  when 
P  Here,  the  FT  estimator  indeed  performs  very 
similar  to  the  ML  estimator,  and  its  MSE  approaches 
the  CRB  together  with  that  of  the  ML  estimator.  In 
addition,  this  case  exemplifies  a  scenario  in  which 
only  one  threshold  effect  is  observed  and  the 
ML  estimator  never  achieves  the  standard  quantum 
limit,  because  at  Nj  =  Nf  the  CRB  has  already 
deviated  from  its  standard-quantum-limited  asymp¬ 
tote.  It  is  evident  from  these  examples  that  determin¬ 
ing  and  Nu  is  critical  to  identifying  the  regime  in 
which  the  ML  estimator  performs  best. 

A.  Lower  Threshold  Nf 

We  will  restrict  our  derivation  to  the  case  of  linear 
frequency  modulation.  Our  approach — motivated 
by  prior  literature — ^will  be  to  first  introduce  an  ana¬ 
lytic  model  that  approximates  the  MSE  of  the  ML 
estimator,  and  then  to  derive  Nf  from  this  model 
[22-24] .  State-of-the  art  photon-counting  detectors 
yield  typically  better  than  kilohertz-rate  dark 
counts,  which  implies  no  more  than  low  tens  of  dark 
counts  per  realistic  integration  windows 
(~1-10  msec).  Because  we  observe  that  the  threshold 
often  occurs  at  much  greater  photon  numbers,  we 
assume  N^  «  0  in  subsequent  derivation.  We  also  as¬ 
sume  that  Nf  occurs  in  a  regime  where  the  photon 
flux  is  low  enough  to  warrant  the  a  «  0  approxima¬ 
tion.  We  begin  by  assuming  that  all  possible  out¬ 
comes  can  be  grouped  into  two  classes  of  events: 
(1)  typical  events,  i.e.,  those  outcomes  that  are  suffi¬ 
ciently  informative  and  yield  average  performance 
close  to  the  CRB,  and  (2)  atypical  events,  i.e.,  those 
outcomes  that  give  little  to  no  information  and  yield 
average  performance  close  to  a  random  uniform 
guess  of  to  £  T.  Based  on  this  classification  we  write 
the  MSE  model  for  the  ML  estimator  as 


is  the  MSE  of  a  uniform  random  guess.  Here  we  have 
used  Ti  =  miuj^gT-to  and  T'2  =  aiax^^gT-to  to  denote  the 
end  points  of  T.  From  the  Eq.  (2^  model  we  define 
Nf  as  the  Nj  value  that  satisfies 


=  (l-P^)CRB.  (25) 

Thus,  good  agreement  between  the  model  and  the  si¬ 
mulated  performance  hinges  on  an  accurate  charac¬ 
terization  of  P^.  When  the  total  photoelectron  count 
is  zero  or  one,  the  ML  estimate  is  uniform  over  T,  cor¬ 
responding  to  atypical  events.  So,  we  can  write 


Pa  =  Pi^  <2)  +  P{n  >  2)P^|2,  (26) 

where  n  =  ni  +  0,2  is  the  total  number  of  counts  from 
both  detectors  within  the  observation  window, 
P{n  <  2)  =  (1  -I-  r]Nj)e~'Ni^  and  P^|2  denotes  the 
atypical  event  probability  conditioned  on  two  or  more 
arrivals.  Unfortunately,  it  is  impractical  to  ascertain 
the  precise  subset  of  outcomes  that  are  atypical. 
Therefore,  we  resort  to  deriving  a  tractable  approxi¬ 
mation  to  P^|2  in  the  ^  1  and  Nj  »  1  limit,  and 
then  using  it  universally. 

To  evaluate  P^|2,  recall  from  Eq.  (^)  that  for  ^  1 
and  (pit)  =  NftjT  the  cost  function  has  an  effective 
resolution  width  equal  to  Aik! Nf.  Dividing  T  into 
^bin  +  1  bins  with  centers  at  t^  =  Ti+  kAnj Nf,  for 
k  =  0,  ...,nbin)  we  approximate  P_4|2  by 

Pa\2  =  l-p(f]  Cih^)  >  C(fy)) ,  (27) 

where  C{t)  is  the  cost  function  given  in  the  argument 
of  the  right-hand  side  of  Eq.  (^),  and  k*  denotes  the 
bin  index  that  contains  the  true  value  ^o-  In  other 
words,  a  typical  event  occurs  if  the  cost  function  is 
peaked  at  the  bin  encompassing  the  true  ^o-  A  trac¬ 
table  expression  for  the  right-hand  side  of  this 
equation  follows  from  assuming  C(4)  are  jointly 
Gaussian  random  variables.  The  mean  of  C(4)  is 
(C(4))  =  {pNiPi^ /2)5k  k,,  and  the  covariance  is  given 
by  cov(C(4),C(fy))  =  {riNiP‘^/2)5kfi  where  <5  is  the 
Kronecker  delta  function.  Consequently,  via 
Eq.  (27)  we  obtain 


P_4|2  »  1  ■ 


■Q[V^iP/2) 


^bin 


(28) 


((^0  -  ^o)">  =  cj\Pa  +  CRB(1  -  P^),  (23) 

where  the  CRB  is  given  by  Eq.  (A6),  but  with  N^  =  0 
and  a  =  0  it  simplifies  to  the  Eq.  (6)  expression.  P^ 
denotes  the  probability  that  the  observed  photon- 
arrival  sequence  belongs  to  the  atypical  set,  and 


2  _  (T2  -  Tif 
^  A  19 


-I- 


T2  +  T1 

2 


(24) 


where  Q(;c)  =  /“ j \/^.  Substituting  Eq.  (28) 
back  into  Eq.  (26)  yields  Pa- 
Having  evaluated  all  of  the  components  in  Eq.  (25), 
we  can  solve  for  N f.  We  perform  several  simplifying 
approximations  to  derive  a  tractable  analytic  expres¬ 
sion.  First,  via  Eq.  (6)  we  use  CRB"^  »  rfNiP'^e^. 
Next,  we  use  Pa  «  Pa\2  in  Efi-  (26),  and  (1  -  x)”  «  1  - 
nx  and  Q(x)  «  e"^^/^/2  in  Eq.  (^),  to  solve  for  the 
breakdown  threshold  as 
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where  W_-^{y)  is  the  Lambert-W  function  along  the  -1 
branch,  formally  defined  as  the  inverse  of  the  func¬ 
tion  y  =  xe^  for  x  <  -1  [25], 

In  Figs.  3  and  4  we  have  plotted  the  N f  predicted 
by  Eq.  (^),  as  well  as  the  MSE  model  from  Eq.  (23) 
and  the  simulation-based  MSE.  In  all  three  cases, 
the  approximation  to  N f  matches  very  well  with 
the  knee  in  the  MSE  model,  indicating  that  it  is  a 
good  analytic  approximation  to  the  transition  in 
the  model.  In  comparing  the  model  to  the  simulated 
MSE  results,  we  find  that  'N f  is  higher  than  the 
breakdown  threshold  observed  in  simulation,  but 
within  a  factor  of  three.  Therefore,  Eq.  (^)  offers  a 
tractable  approximation  to  the  breakdown  threshold 
that  conservatively  predicts  the  Ni  >  Nf  region  in 
which  the  ML  estimator  may  approach  the  standard 
quantum  limit. 

B.  Upper  Threshold  Nu 

Dead  time  impairs  the  estimation  performance  when 
the  average  arrival  rate  of  the  incident  photons  sig¬ 
nificantly  exceeds  1/t,  such  that  the  detectors  cannot 
respond  to  many  information-carrying  incident  pho¬ 
tons.  As  a  result,  the  CRB — hence  the  MSE  of  the  ML 
estimator — starts  to  saturate  when  A'/  >  A„.  To  find 
the  threshold  at  which  this  occurs,  we  return  to  the 
CRB  asymptotes  stated  in  Eq.  (IR)  and  evaluate 
the  Ni  at  which  the  latter  two  asymptotes  intersect. 
The  precise  photon-flux  threshold  at  which  the  CRB 
begins  to  deviate  from  the  standard  quantum  limit  is 
heuristic  in  nature;  however,  motivated  by  the  Fig.  2 
plot  we  define  it  to  be  a  decade  prior  this  intersection 
value,  i.e.. 


_  to)T 
~  5K2ifi,to)T’ 


(30) 


where  K2  and  are  given  in  Eqs.  (Ifi)  and  (IR),  re¬ 
spectively.  This  expression  can  be  evaluated  numeri¬ 
cally  for  a  given  and  but  to  obtain  a  simple 
approximation  we  use  «  A2,  such  that 


Here  the  maximum  is  due  to  the  fact  that  if  the  first 
term  inside  the  brackets  is  less  than  Nf,  then  atypi¬ 
cal  events  dominate  the  performance  and  the 
dead-time-induced  blocking  behavior  is  not  visible 
until  N[  >  Nf. 

The  Nu  we  have  defined  above  is  plotted  in  the  two 
simulations  presented  earlier  in  Fig.  4.  In  both  cases 
it  predicts  the  regime  in  which  dead  time  deterio¬ 
rates  the  MSE  of  the  ML  estimator  reasonably  well. 
In  Fig.  4(a)  the  three  performance  regimes  are  dis¬ 
tinctly  identifiable,  and  the  MSE  approaches  the 
standard  quantum  limit  only  when  Nf  <  Ni  <  Nu, 


as  expected.  In  Fig.  4(b),  on  the  other  hand,  we  illus¬ 
trate  anNf  =  Nu  case:  when  the  incident  photon  fiux 
is  sufficiently  high  to  emerge  from  the  regime  in 
which  atypical  events  dominate,  it  has  already  be¬ 
come  too  high  for  the  dead  time  of  the  photodetectors. 
Consequently,  the  MSE  transitions  from  the  first 
breakdown  regime  directly  into  the  dead-time- 
blocked  regime,  and  never  achieves  the  standard 
quantum  limit. 

In  summary,  the  ML  estimator  can  achieve  the 
standard  quantum  limit  with  the  optimal  1/A/  scal¬ 
ing  of  the  MSE,  provided  that  the  total  mean  photon 
number  incident  on  the  two  detectors  satisfies 
Nf<Ni  <  Nu,  where  A,f  is  given  in  Eq.  (^)  and 
Nu  is  given  in  Eq.  (^).  If  we  have  A„  =  Nf,  then 
the  standard  quantum  limit  is  not  attainable. 

4.  Experimental  Results 

We  have  performed  a  proof-of-concept  FMCW  rang¬ 
ing  experiment  using  linear  frequency  modulation 
and  photon-counting  Geiger-mode  avalanche  photo¬ 
diodes  (APDs).  Here  we  report  on  the  performance 
of  the  ML  estimator  and  discuss  insights  gained  from 
this  experiment. 

Figure  5  provides  a  block  diagram  of  the  experi¬ 
mental  setup.  A  linear  sawtooth  frequency  modula¬ 
tion  with  rate  10.42  GHz /ms  and  period  3.6  ms  is 
applied  to  a  1539  nm  center-wavelength  distributed 
feedback  laser,  using  active  phase  locking.  10%  of  the 
output  power  is  used  in  this  feedback  loop.  Further¬ 
more,  only  the  central  linearized  portion  of  each 
period  is  used  in  collecting  the  experimental  data, 
to  avoid  transients  that  occur  around  the  sharp  fre¬ 
quency  transitions.  The  remaining  output  power 
(90%)  passes  through  several  attenuators,  including 
a  variable  attenuator  for  power  adjustment,  and 
99/1  splitters  for  power  monitoring.  An  in-line  polar¬ 
izer  is  used  to  maintain  polarization.  This  source 
field  is  then  split  into  the  signal  and  reference  arms 
via  a  50/50  splitter.  The  signal  field  passes  through  a 
circulator  and  then  is  sent  out  to  the  target  mirror, 
which  is  placed  approximately  0.8  m  away  on  the 
optical  bench.  The  return  field  is  recoupled  into 
the  circulator  and  is  recombined  with  the  reference 
field  on  a  50/50  splitter.  The  resulting  beat  frequency 
is  approximately  12.7  kHz,  indicating  a  relative 
path  length  of  approximately  18  cm.  The  outputs 
of  the  beam  splitter  are  incident  on  two  matched 
Geiger-mode  photon-counting  APDs,  which  have 


90/10  To  Locking  Short  Delay 


Fig.  5.  (Color  online)  Block  diagram  of  FMCW  ranging 
experiment. 
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quantum  efficiency  of  ??  =  0.1,  dark  noise  level 
Arf  «  lOOOcounts/s,  and  dead  time  t  «  1  |is.  The  out¬ 
put  photoelectrons  in  each  detector  are  time  tagged 
via  a  time-to-digital  converter  unit  over  a  period  of 
T  =  3.4  ms.  The  reference-arm  path  includes  a  short 
delay  line  for  fine  adjustments  to  the  beat  frequency, 
and  an  in-line  polarizer  to  maintain  polarization. 

The  experiment  was  run  at  17  flux  levels  spanning 
approximately  from  80  to  24,000  mean  incident  pho¬ 
tons  per  period  [^],  and  at  each  flux  rate  approxi¬ 
mately  8000  independent  data  sets  were  collected. 
Two  estimators  were  applied  to  each  set,  the  ML 
estimator  from  Eq.  (17)  and  the  FT  estimator  from 
Eq.  (^),  but  with  minor  modification  to  account 
for  the  unknown  carrier-induced  phase  shift  in 
the  ML  estimator  xj/Q  and  to  were  estimated  jointly, 
whereas  in  the  FT  estimator  the  real  part  of  the 
FT  was  replaced  with  its  magnitude  to  render  the 
estimate  independent  of  i//o- 

At  each  photon-flux  level  the  square  error  in  the 
estimates  was  averaged  to  calculate  an  MSE.  A  fair 
assessment  of  the  estimator  performance  requires  an 
accurate  knowledge  of  the  true  range  to  the  target. 
A  calibration  run  with  high  photon  flux  (yielding  ap¬ 
proximately  1600  registered  counts  per  integration 
window  corresponding  to  approximately  24,000  inci¬ 
dent  photons)  was  performed,  and  the  average  of  the 
ML  estimate  and  the  FT  estimate  at  that  photon  flux 
was  assumed  as  the  true  range  to  the  target  [27]. 
Figure  6  shows  the  experimental  RMSE  obtained 
with  the  ML  estimator  and  the  FT  estimator,  as  a 
function  of  rjNi.  The  CRB  and  the  analytic  error 
model  for  the  ML  estimator  are  also  plotted  for  com¬ 
parison.  Both  estimators  demonstrate  good  agree¬ 
ment  with  the  first  threshold  rjNf  predicted  by  the 
theoretical  model.  In  addition,  as  predicted  by  our 
analysis  and  simulations  in  Section  3,  the  ML 
estimator’s  MSE  converges  to  the  CRB  for 
Nj  >  Nf,  whereas  the  FT  estimator’s  MSE  is  approx¬ 
imately  twice  the  CRB  (the  observed  ratios  of  the  FT 
estimator’s  MSE  to  that  of  the  ML  estimator  for 


Fig.  6.  (Color  online)  Performance  of  the  ML  estimator  and  the 
FT  estimator,  plotted  as  a  function  of  ijNj.  The  ML  error  model, 
the  CRB  (dashed  red),  and  the  asymptotic  RMSE  of  the  FT  esti¬ 
mator  (dash-dotted  blue)  are  plotted  for  comparison. 


Ni  >  Nf  ranges  from  1.39  to  2.56).  The  divergence 
of  the  CRB  from  the  standard  quantum  limit  above 
t]Nu  is  not  too  pronounced  for  the  range  of  photon 
fiuxes  plotted  in  the  figure,  but  we  observe  that 
the  ML  estimator’s  MSE  continues  to  demonstrate 
good  agreement  with  the  CRB  above  this  threshold 
as  the  performance  transitions  to  the  dead-time- 
dominated  regime.  There  are  two  data  points  close 
to  the  Nf  threshold  in  which  the  ML  estimator 
appears  to  achieve  slightly  better  performance  than 
the  CRB,  which  we  conjecture  is  a  cumulative  effect 
of  small  errors  in  estimates  of  the  chirp  rate,  the  ob¬ 
servation  period,  the  mean  incident  photon  fiux,  and 
the  residual  statistical  error  in  estimating  an  ensem¬ 
ble  average  from  a  finite  sample  set. 

5.  Discussion 

In  this  paper  we  have  provided  detailed  analysis 
of  optimal  range  estimation  using  frequency- 
modulated  ideal  laser  light  and  continuous-time 
photon-counting  detectors  that  output  time-tagged 
photon  arrivals.  Our  analysis  has  centered  around 
three  themes:  understanding  the  best  possible  per¬ 
formance  via  a  study  of  the  CRB,  characterizing 
the  performance  of  the  ML  estimator,  and  experi¬ 
mentally  validating  the  theory  with  a  proof-of- 
concept  experiment.  The  theory  that  we  have 
developed  is  sufficiently  broad  to  provide  a  unified 
treatment  of  some  of  the  important  recent  advances 
in  employing  different  frequency  modulations,  as 
well  as  varying  signal  and  reference  field  strengths. 
It  can  also  be  applied  to  other  optical  remote-sensing 
instruments  where  frequency  estimates  are  calcu¬ 
lated  from  coherently  detected  fields,  e.g.,  coherent 
Doppler  LIDAR.  Our  CRB  results  reveal  that  the 
best-case  MSE  scales  as  1/A^/.  When  the  reference- 
arm  mean-photon  number  is  comparable  to  that  of 
the  target  return,  the  primary  feature  of  the  (base¬ 
band)  frequency  modulation  that  impacts  the  CRB 
is  its  average  energy,  which  has  permitted  us  to 
quantify  the  gains  from  employing  different  modula¬ 
tions.  When  photodetector  dark  counts  and  dead 
time  are  included,  the  CRB  admits  to  three  distinct 
regimes:  (1)  if  dark  counts  dominate  over  signal 
counts,  the  CRB  scales  as  N^/Nj;  (2)  if  the  signal 
counts  are  high  enough  to  dominate  over  dark  counts 
but  are  not  so  high  that  dead  time  blocks  the  detec¬ 
tors,  the  standard  quantum  limit  prevails;  and  (3)  if 
dead  time  blocks  observation  of  most  incident  photon 
arrivals,  the  CRB  is  constant  with  W/.  Simulating  the 
ML  estimator  has  revealed  two  thresholds,  Nf  and 
Nu,  such  that  when  Nf  <  Ni  <  Nil,  the  performance 
is  approximately  equal  to  the  standard  quantum 
limit,  but  otherwise  the  performance  is  bounded 
away  from  this  optimal  scaling.  We  have  derived  ana¬ 
lytic  approximations  to  these  thresholds  when  linear 
frequency  modulation  is  used,  which  identifies  the 
optimal  performance  regime  in  terms  of  critical  sys¬ 
tem  parameters.  Finally,  we  have  performed  an  ex¬ 
perimental  proof  of  concept  with  equal  reference 
and  target-return  field  strengths,  demonstrating 
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that  the  ML  estimator  behaves  as  predicted  by  our 
analysis,  and  that  it  mildly  outperforms  the  RMSE 
of  the  FT  estimator  (approximately  a  factor  of  a/2) 
at  the  standard  quantum  limit. 

While  the  ML  estimator  has  desirable  optimality 
properties,  it  also  has  some  shortcomings  that  should 
be  addressed.  First,  it  often  has  higher  implementa¬ 
tion  complexity,  as  the  cost  function  to  be  maximized 
may  not  be  efficiently  computable.  In  our  simulations 
the  ML  estimator  required  seconds  to  generate  an 
estimate,  while  the  FT  estimator’s  evaluation  was 
subsecond  (we  did  not  attempt  to  optimize  either 
algorithm,  but  the  FT  estimator  benefited  from  the 
computationally  efficient  fast  FT).  Second,  the  ML 
estimator  is  more  sensitive  to  auxiliary  parameters 
than  the  FT  estimator,  such  as  p  and  ij/Q.  Simplifica¬ 
tions  to  the  ML  estimator  are  often  sought  to 
speed  up  computation  time,  reduce  complexity,  or 
improve  robustness  at  the  expense  of  some  loss 
in  performance.  Although  here  we  have  not  ad¬ 
dressed  these  trades,  the  FT  estimator  does  emerge 
as  a  very  good  candidate  for  linear  frequency 
modulation. 

In  a  field-deployed  ranging  system,  several  addi¬ 
tional  uncertainties  resulting  from  the  environment 
are  significant.  In  particular,  in  this  paper  we  have 
not  considered  the  effects  of  rough-surface  scattering 
off  of  the  target,  atmospheric  turbulence-induced 
amplitude  and  phase  fluctuations  in  the  target- 
return  signal,  ambient  background  radiation,  photo¬ 
detector  timing  jitter,  and  calibration  errors. 
The  theory  presented  in  this  paper  could  be,  in  prin¬ 
ciple,  extended  to  encompass  these  effects  by  deriv¬ 
ing  the  proper  joint  probability  density  function  of 
the  photon-arrival  times.  However,  this  may  prove 
too  complex  for  tractable  and  insightful  analytic 
results. 

In  summary,  we  have  provided  a  rigorous  and  de¬ 
tailed  analysis  of  FMCW  ranging  with  ideal  laser 
fields  and  photon-counting  detectors.  We  have 
derived  the  CRB,  including  the  impact  of  dark  counts 
and  dead  time.  We  have  shown  that  the  ML  estima¬ 
tor  can  attain  the  optimal  CRB  scaling  of  1/Nj  if 
the  mean  received  photon  number  satisfies 
N,<Nj  <  Nu-  We  have  verified  the  performance 
of  the  ML  estimator  via  a  proof-of-concept  experi¬ 
ment  and  shown  that  it  performs  better  than  the  con¬ 
ventional  FT  estimator  as  predicted  by  analysis. 


Appendix  A:  Derivation  of  Cramer-Rao  Bound  and  ML 
Estimator 

We  derive  the  CRB  and  the  ML  estimator  for  general 
dark  count  rate  and  dead  time.  Setting  these  param¬ 
eters  to  zero  yields  the  near-ideal  results  in  Section  2. 
Consider  a  pair  of  matched  photon-counting  detec¬ 
tors  in  the  Fig.  1  setup,  having  quantum  efficiency 
t],  mean  dark  rate  (in  photoelectrons/seconds) 
and  dead  time  t  (in  seconds).  In  the  following  deriva¬ 
tion  we  will  use  =  A^rf(^;^o)  to  keep  our 

notation  compact. 


Taking  the  logarithm  in  Fq.  (11)  yields 


l0gp(ni,Ui,  ...,Un^,n2,Vi,  ....Vn^'to) 

=  log  +  log  f2  +  Cl,  (Al) 


where  we  have  omitted  a  constant  term  independent 
of  ^0-  In  the  limit  that  the  rate  functions 
vary  slowly  over  the  duration  t,  i.e.,  when 
|dD(^;^o)/d^|T  ^  1,  we  have 

/TI2 
T/2 

Substituting  Eq.  (A2)  into  Eq.  (Al),  and  that  into 
Eq.  (5),  yields 


CRB 


T/2 

( 

T/2 


^mAt)  ^mAt) 


-  A,„At) 


(A3) 


where  we  have  used  the  shorthand  notation 
j^mAC  =  ^^mAt)/^to  and  A^,d(C  =  d?A^Ai)/^tl- 
The  last  step  of  evaluating  the  CRB  is  to  find 
In  general,  this  expectation  is  difficult  to  de¬ 
rive  for  self-exciting  stochastic  processes.  However, 
when  the  rate  function  varies  slowly  relative  to 
the  dead-time  interval,  it  has  been  shown  that  the 
average  photocurrent  is  expressible  as  [28] 


{imit)) 


^mAt) 
1  + 


(A4) 


This  mean  arrival  rate  has  intuitive  asymptotes. 
When  «  1,  we  have  «  A^A*)’  indi¬ 

cating  that  the  impact  of  dead  time  is  negligible. 
On  the  other  hand,  if  TA^^(t)  »  1,  the  detector  is 
continually  blocked  by  arrivals  occurring  every  t 
seconds,  so  the  mean  arrival  rate  converges 
to  {iAt))  «  1/r. 

Substituting  Eq.  (A4)  into  (A3),  recognizing  that 
Ai(0  =  -A2(0  and  Ai(0  =  -A2(t),  and  some 
straightforward  algebraic  manipulation  results  in 


CRB-i=  P^dtAL(0[Vv - , .  ^  ^  . 

J-T/2  Am,(i(0(l  +  TAm,d(0). 

(A5) 


Finally,  substituting  the  Eq.  (10)  rate  functions  into 
this  expression,  we  arrive  at 
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CRB 


I 


■T/2  plsinHD(f,to)) 


T/2  i 


■  l^icos^(D(t-,to)) 


1  +  a[l  +  cos^(I)(l;  ^o))] 


1  +  2a  +  a^[l  -  cos^(I)(t,  ^o))]. 


^o) 

(A6) 


where  a  =  r/Nj  dT/{2T).  The  Eq.  (13)  asymptotes  fol¬ 
low  from  considering  the  corresponding  limits  of  this 
expression.  To  obtain  the  standard  quantum  limit 
stated  in  Eq.  (6),  we  set  =  0  and  a  =  0. 

To  derive  the  ML  estimator,  we  start  with  the  log¬ 
arithm  of  the  probability  density  function  that  we 
have  derived  above.  Substituting  Eqs.  (Al)  and 
(A2)  into  Eq.  (17),  followed  by  some  straightforward 
simplifications  of  the  expressions,  yields  the  linear 
filtering  form  of  Eq.  (18)  with  the  impulse  responses 
specified  in  Eq.  (19). 
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